Adiabatic pumping through an interacting quantum dot with spin-orbit couphng 
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We study adiabatic pumping through a two-level quantum dot with spin-orbit coupling. Using 
a diagrammatic real-time approach, we calculate both the pumped charge and spin for a periodic 
variation of the dot's energy levels in the limit of weak tunnel coupling. Thereby, we compare the 
two limits of vanishing and infinitely large charging energy on the quantum dot. We discuss the 
dependence of the pumped charge and pumped spin on gate voltages, the symmetry in the tunnel- 
matrix elements and spin-orbit coupling strength. We identify the possibility to generate pure spin 
currents in the absence of charge currents. 
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I. INTRODUCTION 

A central issue in the field of spintronics is the design 
of spin-based electronic devices. ^'^ They may involve fer- 
romagnets or external magnetic fields to control the spin 
degree of freedom.^ But recently, also all-electric spin- 
tronic devices have gained interest.'^ They rely on spin- 
orbit (SO) interaction, the strength of which is tunable 
via external gates in semiconductor heterostructures,^'^ a 
basic requirement for the realization of a spin field-effect 
transistor.'^ ""^^ 

A spin-polarized current in a semiconductor can be 
generated by spin injection. ^^^^ Here, we focus on an 
alternative route that relies on pumping. By varying the 
parameters of a mesoscopic system periodically in time, 
a finite charge or spin current can be sustained. Experi- 
mental studies have investigated charge pumping in sev- 
eral mesoscopic devices. ^^^^ Spin pumping has been ex- 
perimentally realized in the presence of an external mag- 
netic field. Theoretical studies of spin pumping involve 
external magnetic fields, ferromagnetic leads, ^^^^ and 
also SO coupling.^^^^ 

In the present paper, we consider the minimal model 
that contains SO interaction: a quantum dot with two, 
spin-degenerate orbital levels. Such a two-level quantum 
dot with more than two leads has been suggested as a 
spin filter. We focus on the adiabatic limit of pumping, 
i.e., the parameters are varied slowly in time compared 
to the dwell time of the mesoscopic system. Adiabatic 
pumping of charge and spin through such a two-level dot 
has been considered in the limit of vanishing charging 
energy. It was found that this system can act as an 
all-electric spin battery^ i.e., a finite spin current can be 
achieved without ferromagnets by electrically controlling 
the dot parameters. For specific symmetries in the tunnel 
coupling of the dot to the leads even pure spin currents 
have been suggested. From the analysis of Ref. 28, which 
was based on a scattering-matrix approach, ^^^^ it is not 
clear whether and how the conclusions can be transferred 
to quantum dots with non-vanishing Coulomb interac- 
tion. To answer this question is the main goal of the 



present paper. 

In order to take the Coulomb interaction into account, 
we use a diagrammatic real-time approach, ^^^^ that al- 
lows for arbitrary strengths of the Coulomb interaction. 
We focus on the limit of weak tunnel coupling, for which 
we perform a systematic perturbation expansion to low- 
est order. To emphasize the role of Coulomb interac- 
tion, we compare the limit of vanishing Coulomb interac- 
tion with the limit of an infinitely large charging energy. 

The paper is organized as follows. In Sec. II we intro- 
duce the model that describes SO interaction in a two- 
level quantum dot with Coulomb interaction. Sec. Ill 
deals with the technique to calculate the pumped charge 
and pumped spin during one pumping cycle. To study 
the dependence of the pumped charge (spin) on the four 
tunnel-matrix elements in a transparent way, we intro- 
duce in Sec. IV an isospin representation of the orbital 
degree of freedom. Finally, in Sec. V we present the re- 
sults for the pumped charge and pumped spin. 



II. MODEL 

We consider a quantum dot with two spin-degenerate 
orbital levels \acr) (with labels a = 1,2 for the orbital 
and <j =^7^ foi" the spin), tunnel coupled to the left (L) 
and the right (R) lead. The system is described by the 
Hamiltonian 



H = H 



dot 



H 



lead 



(1) 



Here, is the Hamiltonian of the isolated dot, i^^^^^ 

of the leads, and i^*^" of the tunneling between dot and 
leads. 

The Hamiltonian for the isolated quantum dot con- 
tains two parts. The single-particle contribution for the 
two orbit als a with energy 6q;, which are coupled by SO 
interaction, can be cast in the 4x4 matrix 



(2) 
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FIG. 1. (Color online) Energy scheme of the two- level quan- 
tum dot. The two orbital, spin- degenerate levels can be varied 
in time. They are tunnel coupled to the left (L) and the right 
(R) lead, with tunnel-matrix elements Vx^- The leads have 
the same chemical potential fi. 

for the basis {|1 1) , |1 1) , |2 t) , |2 |)}, where the spin 
quantization axis is chosen arbitrarily. Here, a denotes 
the vector of Pauli matrices, ctq is the identity matrix, 
and OigQ IS a real vector describing the SO coupling. The 
matrix in Eq. (2) has the most general form that al- 
lows time-reversal symmetry. It has been used in the 
context of pumping, and was also recently applied to 
electron-transport in the presence of a magnetic field^^ 
and to study the Josephson current through a double-dot 
structure. In the following, we choose the spin quanti- 
zation axis parallel to ago, so that the matrix becomes 
diagonal in spin space. 

The second part of the dot Hamiltonian accounts for 
the charging energy Ec{N — n^)^, where N is the total 
number of dot electrons and Ug an external gate charge. 
Without loss of generality, we can choose Ug = 1/2 (any 
other value can be achieved by a constant shift of the 
energies Ca)- This leads (up to an additive constant) to 
the dot Hamiltonian 

H""' = ^ e„ di,d,, + J2 i^^so (did,, - h.c.) 

aa a 

+ ^ ^ riat'Tlal + ^ ^la^2a' , (3) 

where the operator (ij^^ creates an electron in state \aa) 
and the corresponding number operator is Uaa = d\^a^aa- 
We used the notation a = ±1 for spin parallel (antipar- 
allel) to ago, Q^so = |c^so|, and U = 2Ec. 

The leads are modeled as reservoirs of noninteracting 
electrons, 

F'-d = ^e,cL,c,,,, (4) 

akX 

where c\^^ is the creation operator for an electron with 
spin a and momentum k in lead A. Tunneling between 



dot and leads is described by the Hamiltonian 

ff*""= E ^AacLArf„.+ h.c, (5) 

aakX 

with (spin-independent) tunnel-matrix elements Vxa for 
tunneling between lead A and orbital a. 

Pumping is achieved by varying system parameters pe- 
riodically in time. In this paper, we assume that the en- 
ergy levels Ca {t) can be changed in time via external gates 
capacitively coupled to the system. Although, in princi- 
ple, the external gates may affect also the SO coupling, 
the tunnel couplings and the electro-chemical potential 
of the leads (via parasitic capacitances), we assume all 
these parameters to be constant in time. 

We focus on the regime of adiabatic pumping, which 
is achieved for pumping frequencies Q smaller than the 
inverse of the dwell time. This is valid for ^7 <C T, where 
r is the tunnel-coupling strength, F = X^Aa^^aa, with 
^Xaa' = pVxa'Vxa' The density of states p is assumed 
to be flat and equal for the left and right lead. We choose 
a gauge where all four tunnel-matrix elements are real. 

To study the effect of Coulomb interaction, we com- 
pare results for the limit of noninteracting (/7 = 0) and 
infinitely strong interacting (U = oo) electrons on the 
dot. In the latter case, the total number of electrons in 
the quantum dot can only be zero or one. 



III. METHOD 

To calculate the pumped charge and pumped spin, 
we use a diagrammatic real-time approach to adia- 
batic pumping through quantum-dot systems. For the 
present context, we extend the analysis of Ref. 39 to allow 
for a time-dependent transformation of the basis states. 
This is necessary since the SO coupling couples time- 
dependent orbital levels, which, in turn, makes the dot 
eigenstates time dependent. 

We start in Sec. HI A with the kinetic equation for 
the reduced density matrix in its general form, which de- 
scribes the time evolution of the dot's degrees of freedom. 
Subsequently, we perform both an adiabatic expansion, 
i.e., a perturbation expansion in the pumping frequency 
(Sec. IIIB) and a perturbation expansion in the tunnel- 
coupling strength (Sec. IIIC) to describe the limit of 
weak tunnel coupling. The pumped charge and pumped 
spin currents to lowest order in F and Vt are derived in 
Sec. HID. Finally, in Sec. HIE, we perform the limit 
of weak pumping which assumes small amplitudes of the 
pumping parameters. 



A. Kinetic equation 

The main idea of the diagrammatic real-time technique 
is based on the fact that the leads are described as large 
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reservoirs of noninteracting electrons which can be inte- 
grated out in order to arrive at a reduced density ma- 
trix p for the dot degrees of freedom only. For a matrix 
representation with matrix elements p^^ = (xi | p' 



.dot 



X2) 



(for the diagonal elements we introduce the notation 
= P^)' convenient to use the eigenstates \xi) with 
corresponding eigenenergies E^. as a basis. For this, we 
employ a time-dependent unitary transformation T act- 
ing on the dot Hamiltonian such that t'^H'^^^T is 
diagonal. 

The time evolution of the reduced density matrix is 
given by the kinetic equation 



dt'W{t,t') p{t') 



(6) 



The bold face indicates tensor notation. The reduced 
density matrix p and T^T are tensors of rank 2, while 
AS and W are tensors of rank 4, i.e., 



{W{t,t')p{t% 



K:^:it^t')pliit')- (7) 



i.e., it collects all contributions to first order in the pump- 
ing frequency Q. The difference in the eigenenergies of 
the isolated dot, AS{t)^ is of instantaneous order, while 
T^T{t) belongs to the adiabatic correction. 

Since both W*^*^ and W'^''^ depend only on the dif- 
ference t — t\ it is convenient to perform the Laplace 
transform F{z) = df e-<^-^'^F{t - t'). Using 

= yV^t^''\z = 0+) and 
, the kinetic equation reads 



the short notations W^*^'^^ 



=0+ 



in instantaneous order, and 



(8) 



-P[^ = (W[^--AS ]p[^ 



T^t,p\ 



(i) 



(9) 



Xi,X2 



for the adiabatic correction. The normalization condition 
for the density matrix is expressed as Trp[*^ = 1 and 

Trp(°) =0. 



The kernel element W^^^i^ {t^t') describes the transition 

from p^{t') at time t' to P'^^{t) at time t. It is given 
by the sum over all irreducible blocks on the Keldysh 
contour which correspond to the described transition. 
The elements of AS are differences of the eigenenergies 



defined as {AS(t)Y^^] 



^X2(^)) ^XlXl^X2X2 



The second term 



T^T{t)^p{t) originates from the time 

dependence of the transformation T, and T denotes the 
time- derivative of T. In the following adiabatic expan- 
sion and the expansion in the tunnel-coupling strength, 
we follow the lines of Ref. 39. 



B. Adiabatic expansion 

In the limit of slow variation of the system param- 
eters, such that the duration of one pumping cycle, 
T = 27r/r^, is much larger than the dwell time of an 
electron in the quantum dot, we can perform an adi- 
abatic expansion of Eq. (6), which is equivalent to an 
expansion of all time dependences around the final time 
t and to systematically keep all contributions that con- 
tain one time derivative. For this, we first do a Taylor 
expansion of the reduced density matrix around the fi- 
nite time t, i.e., p(t') p{t) + (t' - t) ^p{t). Then, we 
expand the kernel and the density matrix in the pump- 
ing frequency, i.e., p{t) p^*^ + Pt^^ and W{t,f) 
W['\t - t') + yV^^\t - t'). The instantaneous order, 
indicated by the index (i), describes the limit where all 
system parameters are frozen at time t. The adiabatic 
correction, labeled by (a), contains one time derivative. 



C. Expansion in the tunnel-coupling strength 

In addition to the adiabatic expansion, we perform a 
perturbation expansion in the tunnel-coupling strength 
F. For a systematic expansion of the kinetic equations, 
we need to analyze the term AS p^*/^) . It vanishes for all 
diagonal matrix elements of . The off-diagonal ma- 
trix elements, associated with coherent superpositions, 
are only nonzero when the superposition is not forbidden 
by conserved quantum numbers and when the energy dif- 
ference of the corresponding states is smaller or of the 
order of F. Therefore, we count all contributing matrix 
elements of AS to be of the order of F. 

The expansion of the kernels, W^*^^^ = 

VV^*^^'"'\ starts to first order in F and the 
instantaneous order of the reduced density matrix to 

zeroth order, p^f^ = X^^o^t*'"^^- properly match 
the powers of F in Eq. (9), the adiabatic correction of 
the reduced density matrix, p^^^ = ^'^^_ip[^'^\ has to 
start to minus first order. 

The instantaneous part of the kinetic equation starts 
to first order in F, 



with normalzation Trp^*'^^ = 1. For the adiabatic cor- 
rection, the expansion of Eq. (9) to lowest (zeroth) order 
in F yields 

Ap(»-°) = (wP4A£)p(--), (11) 



(10) 
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with Trpi' 



(a,0) _ 



0. An other terms appearing on the 6ea{t): 



right-hand side of Eq. (9) are of higher order in F and 
drop out. This is immediately obvious for the last two 



terms in Eq. (9). But also T^T.pl 



Am 



drops out in the 



absence of any bias voltage. In this case, Pt~' is given 
by the equilibrium distribution, which is diagonal with 
matrix elements being determined by Boltzmann factors. 
Since energy differences, A£^, of states for which coher- 
ent superpositions are allowed are of the order of F, the 
difference of the corresponding occupation probabilities 
for these states is also of the order of F and, therefore, 
vanishes in the perturbation expansion. This means that 
the matrix elements 



Xi 



(12) 



vanish for all combinations of xi and X2 which are needed 
in the kinetic equation. 



D. Pumped charge and pumped spin 

The pumped current and the pumped spin current 
from the dot into the left lead are given by 



hit) = e f dt' 

J — oo 



Tr 



dt' Tr 



WQ(i,i')p(t') 
W^(t,t')p(i') 



(13) 
(14) 



respectively. Here, we have introduced WQ/5(t,t') = 
Y.q{Qt ± where W^"^^"^^ {t,t') only con- 

tains those diagrams of W(t,t') in which Qo- electrons 
with spin a enter the left lead, i.e., in which the number 
of lines for lead L and spin a going from the upper to the 
lower contour minus the number of lines from the lower 
to the upper contour is q^. 

Analoguously to the expansion of the kinetic equation, 
we perform the adiabatic expansion and the perturba- 
tion expansion in the tunnel-coupling strength for the 
pumped charge and spin current. To lowest order we get 



j[°'°Ht) =eTr 



it) 



Tr 



Pt 



(a-l) 



Pi 



(a-l) 



(15) 
(16) 



The pumped charge and the pumped spin per pumping 
cycle is obtained by integration, Q = jj" dt I^'^\t) and 
S = dt I^'^\t). The diagrammatic rules to calculate 



W^^'^^ can be found in App. A.^^ ^^'^^'^^ 



E. Weak pumping 



We split the energy of the orbital levels into the time- 
averaged part, Ca = ^ lo' ^a{t)^ and the deviation 



ei(t) 
62 (t) 



ei 

£2 



5e2{t). 



(17) 
(18) 



In the limit of weak pumping, the time-dependent part of 
the pumping parameters is small compared to other typi- 
cal energy scales of the system such as temperature, /c^T, 
and tunnel-coupling strength, F, hence we can expand 
the pumped charge and pumped spin to lowest (bilin- 
ear) order in ei{t) and 62 (t). For adiabatic pumping, the 
pumped charge (spin) is proportional to the area enclosed 
by the path of (ei(t), e2(t)) in the parameter space during 
one pumping cycle. Therefore, a phase difference is nec- 
essary to gain finite pumped charge (spin). The enclosed 
area is given hy rj = dt 6ei{t) dtSe2{t). All results in 
Sec. V are calculated in the weak-pumping limit. 



IV. ISOSPIN TRANSFORMATION 

For each matrix element p^^ that needs to be consid- 
ered (all diagonal ones and those off-diagonal ones that 
describe possible coherent superpositions), there is one 
kinetic equation. It is often convenient to transform the 
reduced density matrix such that only linear combina- 
tions of the p^J's appear, which allows for a straight- 
forward physical interpretation. In the context of spin 
transport through a single-level quantum dot with ferro- 
magnetic leads, it is advantageous to formulate the ki- 
netic equations separately for the occupation probability 
of zero, one or two electrons on the quantum dot and 
the three components of the spin on the dot.^^'^^'^^ The 
vector character of the spin accounts for both a spin im- 
balance along a given axis and the coherent dynamics of 
the accumulated spin. One virtue of such a transforma- 
tion lies in the fact that it is possible to write the kinetic 
equation in a (spin-) coordinate- free form, which does not 
depend on the choice of the spin quantization axis. 

A similar transformation can also be used for the or- 
bital degree of freedom in systems in which coherent su- 
perpositions of the occupation of different orbitals ap- 
pear. These superpositions are conveniently described by 
defining an isospin. This has been done before for sev- 
eral double-dot systems. ^^'^^^^ We will introduce such 
an isospin description now for the system under consid- 
eration. 

In this paper, we focus on the limits of U = and 
= 00. In the first case, = 0, the Hilbert space is 16- 
dimensional, i.e., the reduced density matrix is a 16 x 16 
matrix. However, since we choose the spin-quantization 
axes along the direction of the SO field, the Hamiltonian 
divides into two independent spin channels. As a result, 
the reduced density matrix can be written as a direct 
product of the 4x4 density matrices for spin up and 
spin down, pu^o = (pc/=o)^ {pu=o)i' In the basis 
{|0) , |1) , 1 2) , \d)}a- that corresponds, for each spin a, to 
the occupation of none of the orbitals, of orbital 1, of 
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orbital 2, and of both orbitals, respectively, the reduced 
density matrix reads 



/Po 0- 

P2 
VO PdJ 



(19) 



Note that, in order to keep the notation simple, we put 
the index <j only once at the matrix indicating the a- 
dependence of each of the matrix elements. 

For U = oo the dot is either singly occupied or empty, 
i.e., the Hilbert space is 5-dimensional. The reduced den- 
sity matrix takes the form 




Pu= 











It 

Pit P2t 
2t 

Pit P2t 












pit Pnl 



(20) 



Note that, here, po is the probability that the dot is not 
occupied with either spin, while for /7 = we used (po)cr 
for the probability that the dot is not occupied with spin 
cr, irrespective of the occupation of spin —a. 

To describe the coherent superposition associated with 
the off-diagonal matrix elements, it is convenient to in- 
troduce, for each physical spin, an isospin operator I^j 
with quantum-statistical expectation value I^r = (^a)- 
Choosing the coordinate system for the isospin such 
that I Icr) and |2cr) are the eigenstates of we get 



A I- = i {pI- - pI-) /2, and 



{pia — P2a) /2. Since ultimately we aim at a coordinate- 
free form of the kinetic equations, we abbreviate the z- 
axis chosen above by the normalized vector n, i.e., |lcr) 
and |2cr) are the eigenstates of 1^ • n. 

The isospin direction n characterizes the eigenstates of 
the isolated quantum dot in the absence of SO coupling. 
The SO coupling, however, couples the two orbitals. As 
a consequence, the dot eigenstates 



la) 

2a) 



(21) 



for single occupation with spin cr are linear combinations 
of the two orbitals |lcr) and |2cr), given by the transfor- 
mation 



1 



^ + Ae iaaso 
icraso ^ + Ae 



(22) 



The corresponding eigenenergies are E± = e±^, with the 
mean dot level e = (ei + e2) /2 and ^ = y^Ae^ + a'^^. The 
transformation depends on the spin cr, the level spacing 
of the both orbitals, Ae = (ei — e2)/2, and the strength of 
the SO coupling, a^o- We consider the regime where Ofgo 
and Ae are of order F. Therefore, the level spacing 2^ of 
the eigenenergies E± is also of order F. As we pump on 



FIG. 2. (Color online) Scheme of different relevant isospin 
quantization axes. The vector n represents the quantization 
where the orbital levels \la) and |2cr) are the eigenstates of 
the Iz operator of the isospin. The two axes fia are the quan- 
tization axes where the eigenstates of Iz are the eigenstates 
of for single occupation. In a ferromagnetic analogy, the 
leads are fully isospin polarized along the axes mx. 



both energies, ei{t) and e2(t), the transformation Tcr(t) 
and the eigenenergies E±{t) are time dependent. 

The unitary transformation T^r corresponds to a rota- 
tion about the x-axis with the spin-dependent angle 



(23) 



in isospin space. This means that the dot eigenstates 
|±cr) are eigenstates to the isospin projection la-fia- along 
the direction n^r that is obtained from n by the above 
mentioned rotation (see Fig. 2). 

The tunneling Hamiltonian couples the lead-electron 
states to both orbitals, i.e., to a linear combination of 
I Icr) and |2cr). To diagonalize the tunneling from and to 
lead A, we employ the unitary transformation 



-VX2 Vxi 



(24) 



In isospin space, this transformation corresponds to a 
rotation about the ?/-axis with angle 



2VxiVx2 



(25) 



Applying this rotation on n generates the direction m\ 
(see Fig. 2) which has the following physical interpreta- 
tion: only dot electrons with |+)/^.^^ isospin projection 

along Jcr • ''^A couple to reservoir A, while the \—)f 
isospin projection is decoupled from the lead.^^ There- 
fore, in a ferromagnetic analogy, the leads are full isospin 
polarized with polarization along m\. 

Including this isospin in the formulation of the kinetic 
equation, the system is fully described by the occupa- 
tion probabilities of the dot and the expectation values 



6 



of the isospins. In the hmit of vanishing Coulomb in- 
teraction, the probabihties describing the occupation of 
the dot with spin a are {po^Ps^Pd)a- ^^r empty, po^ single, 
Ps = Pi and double occupation, p^. We expand the 
kinetic equations to first order in the tunnel coupling. As 
described above, we count both a^o and Ae as one order 
in r. We get for /7 = 



tions, then, read 




(26a) 



(26b) 



where / = /(e) is the Fermi function at energy e. As the 
difference of the eigenenergies, 2^, is of order F, we have 
to drop 2^ in terms which are already linear in F. There- 
fore, the Fermi function, /, depends here only on the 
mean level position, e, since every term which includes 
the Fermi function is linear in F. In the equations for the 
probabilities, the isospin projections along the directions 
defined by the leads enter in the weighted average 



-mL 



(27) 



with Tx = ^^Txaa- The isospin projection direction 
given by the SO coupling, on the other hand, gives rise 
to a precession term about the effective field 



2^. 



(28) 



in the equation for the isospin. This effective field is the 
only place where the SO coupling enters the kinetic equa- 
tions. Equations (26a) and (26b) represent both the in- 
stantaneous order and the adiabatic correction of the ki- 
netic equation. For the first case, one needs to set the left- 
hand side to zero and add the index (i, 0) to the isospin 
and the occupation probabilities on the right-hand side. 
(Note that the instantaneous part of the isospin vanishes 

in lowest order in F, 1^'^^ = 0.) For the the second case, 
we need to add the index (i, 0) on the left-hand side and 
(a, —1) on the right-hand side, respectively. 

In the limit of strong Coulomb interaction, = oo, 
the relevant occupation probabilities are {po^Pt^Pi) with 
Pa = Pia +P2cr being the possibility that the dot is occu- 
pied by a single electron with spin a. The kinetic equa- 




-2/ (l-/)/2 (l-/)/2 
/ -(l-/)/2 
/ -(l-/)/2y 



+ ^(1-/) -(^f^) 



'Pq' 
Pt 

(29a) 



dt 



Po 



f 



h \2 
r i-f 
h 2 



m 



• X {B^ + Bu) . 



(29b) 



In addition to the effective field B^j generated by the 
SO coupling, there is another effective field acting on the 
isospin. The latter appears as a consequence of the in- 
terplay between tunneling and Coulomb interaction. It is 
identical to the exchange field acting on the physical spin 
in quantum dots attached to ferromagnetic leads. ^^'^^'^^ 
In the limit of U = oo, it is given by 



Br 



F 

F 

2^ 



m Re J 

/^^cutoff 



In- 



27r 



Re^ 



1 

2^ 



(3e 



(30) 



where t/cutoff appearing in the second line is the high- 
energy cutoff required to make the integral convergent. 
The digamma function is denoted by ^, and we used 
P = l/kBT. 



V. RESULTS 

In this section, we present the results for the adiabati- 
cally pumped charge (spin) in the weak-pumping regime. 
To analyze the effect of Coulomb interaction, we com- 
pare results for noninteracting electrons, /7 = 0, with the 
limit of strong Coulomb interaction. The latter is real- 
ized by setting U = oo m the Hamiltonian and, thereby, 
suppressing occupation of the quantum dot with more 
than one electron. Furthermore, finite Coulomb interac- 
tion influences the amplitude of the exchange field, Bu^ 
via the the high-energy cutoff, /7cutoff- In all calculations, 
we set t/cutoff = We assume weak tunnel cou- 

pling between quantum dot and leads, F <^ ksT, i.e., 
we restrict the calculation to lowest order in the tunnel- 
coupling strength F. If not stated otherwise, the SO- 
coupling strength is a^o = F/10. 

For = 0, the results of this paper can be compared 
with calculations that include higher orders in F. For 
example, Brosco et al. have studied the two- level quan- 
tum dot with SO coupling and vanishing Coulomb inter- 
action in the limit of zero temperature.^^ Calculations 
exactly in F can be done, e.g., with a scattering matrix 
approach, ^^^^ which is equivalent to an approach that 
is based on a formula relating the pumped current to 
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the instantaneous dot Greens functions. The latter is 
in general extendable to finite interaction. 

In this section, we study the dependence of the pumped 
charge (spin) on various parameters: the strength of the 
SO coupling, Q^so, the tunnel coupling to the leads, Vaq,, 
and the time-averaged dot levels, e^. It is convenient to 
parameterize the latter by the time-averaged mean dot 
level, e = (ei + e2)/2, and the averaged spacing of both 
orbital levels, Ae = (ei — e2)/2. 

In the regime under consideration, the temperature ap- 
pears only in two ways. First, since the mean energy level 
only appears in the combination /3e, the temperature pro- 
vides the energy scale on which variation of the mean 
level energy changes the pumped charge (spin). Second, 
the absolute value of the pumped charge and pumped 
spin are proportional to {kBT)~ . Therefore, all plots are 
normalized accordingly. The dependences of the pumped 
charge (spin) on the other parameters are not affected by 
temperature. 

The tunnel coupling of the two dot orbitals to the left 
and the right lead is defined by four real tunnel-matrix 
elements, V\a- If the tunnel- matrix elements are equal 
for the coupling to the left and right lead, Vi^^ = V^iai 
then, for symmetry reasons, there will be no pumping 
transport via variation of the quantum dot's levels. To 
achieve pumping, the left-right symmetry needs to be 
broken by changing either the magnitude or the sign 
of one the tunnel couplings. We find it convenient to 
parametrize the tunnel- matrix elements by angles (px^ 
which have been introduced in the previous section (see 
Eq. (25)). Then, the tunnel-matrix elements are given by 

the relations Vai = ^ cos ^ and Va2 = \J sin 

For (j)x = 7r/2 both orbital levels are coupled symmetri- 
cally to lead A, i.e. Vai = Va2, and for (j)x = — 7r/2 the 
orbitals are coupled antisymmetrically, Vai = — Va2- The 
necessary condition to get a finite pumped charge (spin) 
is 7^ since (j)L = ^i? (even for Fl 7^ Fr) leads auto- 
matically to an effective one parameter pumping without 
any finite pumped charge (spin) in the adiabatic limit. 



A. Charge and spin pumping 

Motivated by the previous discussion, we first focus 
on a tunnel-coupling configuration with Fl = Fr but 
7^ ^i?, where the pumped charge and pumped spin 
are, in general, finite. Both depend on the mean dot-level 
positions, e and Ae, which is shown in Fig. 3. Those or- 
bital energies of the quantum dot can, in principal, be ad- 
justed by capacitively coupled gate votages. Figures 3(a)- 
3(d) illustrate the pumped charge (spin) for /7 = and 
U = 00 and for a tunnel-coupling configuration where 
the coupling to the left lead is symmetric regarding the 
orbitals, = 7r/2, while the coupling to the right lead 
is given by = 7r/4, i.e., FRi/yR2 = cot7r/8. In case 
that orbital 1 is symmetrically (I^li = Vri) and orbital 
2 is antisymmetrically {Vl2 = — ^^2) coupled to the left 



and right lead, the pumped spin is in general finite while 
the pumped charge vanishes for this configuration. In 
Figs. 3(e) and 3(f) the pumped spin is exemplarily cal- 
culated for — = (/>R = 7r/4, which is equivalent to 
Vli Vri -Vl2 cot tt/S Vr2 cot 7r/8. 

Each plot in Fig. 3 shows a maximum and a minimum 
value. For no Coulomb interaction, the maximum value 
is located at e = (relative to the chemical potential 
/i of the leads). In the limit of strong Coulomb inter- 
action, the extrema positions are shifted to values of e, 
whose order of magnitude is given by the temperature. 
The Ae-position of the maximum pumped charge (spin) 
depends on the tunnel coupling to the leads and the SO- 
coupling strength. Increasing a^o also increases the max- 
imum's position with respect to Ae. Furthermore, the 
pumped charge is in general larger for no Coulomb inter- 
action apart from special tunnel-coupling configurations 
discussed in detail in the next section. That is not sur- 
prising since the Coulomb interaction reduces the pos- 
sible transport channels through the dot by suppressing 
occupations of the dot with more than one electron. 

B. Exchange- field interaction 

Both limits U = and U = 00 show different sym- 
metries with respect to Ae — Ae. In the limit /7 = 0, 
the pumped charge (spin) is exactly antisymmetric in Ae. 
The antisymmetry with respect to (e, Ae) (— e, — Ae) 
originates from the particle-hole symmetry. The anti- 
symmetry in Ae alone, on the other hand, is a non-trivial 
result and only valid for the lowest order contribution 
in F. In the limit of strong Coulomb interaction, the 
symmetry in Ae differs from the U = limit. The ex- 
change field Bu, which interacts with the isospin, leads 
to a contribution of the pumped charge (spin) that is not 
antisymmetric in Ae. Therefore, the antisymmetry is in 
general broken. To point out the symmetry character- 
istics, we study the pumped charge (spin) in two differ- 
ent tunnel-coupling configurations (1) and (2), both for 
Fl = Tr, 

(1) : = ^ , = Y , 

(2) : (/>L = -J, </>R = ^, (31) 

which is equivalent to (1): Vi,i/Vi,2 = cot7r/8, 
^Ri/^R2 = l/>/3 in the first case, and (2): Fli/Fl2 = 
— cot7r/10, l^Ri/1^2 = cot7r/8 in the second one. Tunnel 
coupling (1) and (2) show that the exchange field can af- 
fect the pumped charge and the pumped spin differently, 
and the effect, thus, depends on the tunnel-coupling pa- 
rameters. That is accounted for by Fig. 4, where the 
cut through the contour plot (of Fig. 3 but with coupling 
configurations (1) and (2)) for fixed e is shown. The fixed 
value of e is e = in the limit of vanishing Coulomb in- 
teraction and, for comparison, e = ksT in the limit of 
strong Coulomb interaction. 
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FIG. 3. (Color online) Pumped charge (spin) with finite SO coupling, ago = T/IO, in the [7 = and the U — oo limit depending 
on the time- averaged orbital level positions. There are two sets of coupling parameters with Vl = r^: (a)-(d) = 7r/2, and 
0R = vr/4. (e)-(f) (/)L — — 7r/4, and (/)r = 7r/4. The latter coupling leads to vanishing pumped charge. 



For configuration (1), the exchange field leads to a peak 
located at Ae = which has a nearly symmetric behav- 
ior in Ae. The pumped spin, on the other hand, is still 
approximately antisymmetric in Ae. Furthermore, with- 
out Bu^ the pumped charge (spin) is usually smaller for 
U = 00, compared io U = 0, because of the reduced 
number of transport channels through the dot, but the 
exchange field can enhance the pumped charge. There 



are sets of parameters where the charge transport is even 
larger for finite Coulomb interaction than for /7 = 0. 

For tunnel coupling (2), the symmetric part of the 
exchange-field contribution is less important. The 
pumped charge, in this case, is not dominated by a sym- 
metric behavior as we observed for coupling (1). It is 
rather a shift of the point of zero pumped charge to a 
finite value of Ae similar to the pumped spin. 
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FIG. 4. (Color online) Pumped charge (spin) with finite 
SO coupling, aso = T/IO, in the [7 = and the U = oo limit 
depending on the averaged level-spacing of both orbital lev- 
els. In the U = limit we choose e = 0, which is the position 
of the maximum value. In the limit of strong Coulomb inter- 
action, we use e = /cfiT as an approximation to the position 
of the maximum value. The two sets of coupling parameters 
(1) and (2) are the ones given in the text (see Eq. (31)). 
The vertical dotted lines indicate pure spin pumping. In the 
weak-coupling limit, this is only possible for pumping with 
Coulomb interaction. 



Comparing the exchange-field contribution for config- 
urations (1) and (2), the contribution to the pumped 
spin reaches its maximum where the contribution to the 
pumped charge vanishes, and it is approximately half 
of its absolute maximum value where the contribution 
to the pumped charge has its maximum. For large val- 
ues of exchange-field contribution, near its maximum, the 
pumped charge has a dominant symmetric contribution 
while the exchange-field contribution to the pumped spin 
is in general too small to generate a peak at Ae = 0. 



C. Pure spin pumping 



Pure spin pumping is achieved whenever the pumped 
charge vanishes but the pumped spin remains finite. To 
find such points it is helpful that the pumped charge 
and pumped spin behave differently in the presence of 
Coulomb interaction, as discussed in the previous section, 
and that the pumped charge is more sensitive to symme- 
try in the tunnel-matrix elements than the pumped spin. 
This defines the two strategies to obtain pure spin pump- 
ing: to tune either the orbital energy levels of the dot or 
the tunnel-matrix elements. 



1. pure spin pumping by tuning orbital energies 

For fixed tunnel couplings, we try to tune the orbital 
energies such that the pumped charge vanishes but the 
pumped spin remains finite. As discussed above, this is 
easily possible for strong Coulomb interaction, because 
in this case, the value of Ae at which the pumped charge 
changes its sign is shifted away from Ae = due to 
the exchange field Bu. In absence of Coulomb interac- 
tion (and to lowest order in the tunnel coupling strength), 
this does, in general, not work apart from special cou- 
pling configurations, where the pumped charge vanishes 
independently of the orbital energies, as discussed in the 
next section. The reason is that both the pumped charge 
and the pumped spin are, to lowest order in F, exactly 
antisymmetric in Ae, i.e., the pumped charge and spin 
vanish simultaneously. The comparison between the two 
limits is shown in Fig. 4. The points of pure spin pumping 
are indicated by the vertical dotted lines. Another inter- 
esting feature of the finite difference between the zero- 
points for the pumped charge and the pumped spin is 
the possibility to change the sign of the pumped spin, 
while charge is pumped in the same direction. 



2. pure spin pumping by tuning tunnel couplings 

There are cases in which pure spin current is not only 
possible for special, fine-tuned orbital energies but for all 
values of e and Ae. This is illustrated in Fig. 5, which 
shows the maximum absolute value of the pumped charge 
(spin) in the (e, Ae) parameter space, as a function of the 
coupling parameters (j)x for Fl = Fr and for both limits 
U = and U = oo. The plots can be periodically contin- 
ued. The dotted lines represent coupling configurations 
where the pumped charge and the pumped spin are zero. 

Along the middle dotted line, = pumped charge 
and pumped spin vanish due to left-right symmetry as 
mentioned previously. Here, the tunnel-matrix elements 
are equal for the coupling to the left and the right lead, 
^La = yRa- The dotted zero-lines (/)r = i tt for 
zero pumped charge (spin) only exist for lowest order in 
F; higher-order corrections would lead, in general, to a 
finite pumped charge (spin). The latter conclusion can be 
drawn by comparing with calculations for C/ = which 
are exactly in F, e.g., by means of a scattering matrix 
approach, ^^'^^^^ and it is self-evident that even finite 
Coulomb interaction does not change that significantly. 
Along these lines the tunnel-matrix elements are given 
by "l^Li^Ri = —"^2^2- In any case, these dotted lines 
do no mark good candidates for pure spin pumping since, 
there, charge and spin behave similar. 

The situation is different along the dashed lines. The 
middle dashed line, (/)r = — represents a configuration 
where for each orbital the absolute value of the tunnel- 
matrix elements is the same, but one element of all four 
has an opposite sign, i.e. Vi,i = Fri and Vi,2 = — "^^2 
(or equivalently 1^2). Here, we find (to lowest order 
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FIG. 5. (Color online) Pumped charge and pumped spin in the [7 = and the U — oo limit depending on the orbital-coupling 
configuration for fixed Fl = Fr. The illustrated function shows the maximum value of the pumped charge (spin) in the 
(e, Ae) parameter space for a^o — F/10. The dotted lines represent coupling configurations where the pumped charge (spin) is 
zero. Along the dashed lines, for [/ = 0, the pumped charge is always zero while the spin is still finite. For strong Coulomb 
interaction, U — oo^ only the line = — 0r leads to vanishing pumped charge. 



in F) pure spin pumping for both vanishing and strong 
Coulomb interaction. This generalizes the result found in 
Ref. 28 for the [7 = limit to the limit of strong Coulomb 
interaction. The dependence of the pure pumped spin for 
= = 7r/4 on e and Ae in both Coulomb regimes 
is shown in Figs. 3(e) and 3(f). 

The dashed lines = — i tt (equivalent to 
^Li^Ri = ^L2^R2) indicate a further scenario for pure 
spin pumping to lowest order in F in the = limit. 
For higher orders in F, however, the pumped charge be- 
comes finite. It also becomes finite ior U = oo (and 
lowest order in F) as a consequence of the exchange field 
acting on the isospin. 

How important is the symmetry Fl = Fr? To answer 
this question, we calculate the pumped charge and spin 
for Fl = 2Fr, see Fig. 6. As we see, the dependence of 
the pumped charge and spin on (px changes substantially 



for the pumped charge but not so much for the pumped 
spin. In particular, there are no straight lines with pure 
spin pumping anymore. For /7 = (and to lowest order 
in F), pure spin pumping is still possible on curved lines 
in the parameter space, but not for /7 = oo. There- 
fore, Fl = Fr is a necessary requirement for pure spin 
pumping. 



D. Spin-orbit coupling strength 

The dependence of the pumped charge and pumped 
spin on the SO-coupling strength is visualized in Fig. 7. 
Here, the different functions again show the maximum 
value of the absolute pumped charge (spin) in the (e, Ae) 
parameter space. The upper plot illustrates the pumped 
charge. In general, the SO coupling reduces the pumped 
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FIG. 6. (Color online) Pumped charge and pumped spin in the [7 = and the U — oo limit depending on the orbital-coupling 
configuration for fixed Fl = 2Fr. The illustrated functions show the maximum value of the pumped charge (spin) in the (e, Ae) 
parameter space for a^o — F/10. 



charge during one pumping cycle. It will also decrease 
for large Ae. In both cases, the difference of the eigenen- 
ergies is large so that the pumping mechanism is an ef- 
fective one-level pumping, which does not lead to any 
pumped charge (spin) in the adiabatic regime. For small 
values of a^o compared to F, the effect of the exchange 
field is enhanced. In this regime, the Coulomb interac- 
tion increases the pumped charge. For larger values of 
Qfso, the influence of the exchange field is reduced. The 
pumped spin, on the other hand, is not crucially affected 
by the exchange field. 



The pumped spin, in contrast to the pumped charge, 
vanishes for ago = 0- Therefore, there is an optimal value 
of Qfgo that maximizes the pumped spin (see Fig. 7). This 
value is smaller than F and it depends on the tunnel 
coupling. 



VI. CONCLUSION 



We analyze the possibility to build an all-electric spin 
battery and to generate a pure spin current with a two- 
level quantum dot in the presence of Coulomb interac- 
tion. In the limit of vanishing Coulomb interaction, both 
is possible, as has been demonstrated in Ref. 28. Here, 
we show that this is also possible for the experimentally 
relevant case of a quantum dot with large Coulomb inter- 
action. The Coulomb interaction changes the pumping 
characteristics substantially. In particular, symmetries 
with respect to the orbital energies change as a conse- 
quence of an effective exchange field acting on an isospin 
defined by the orbital level index. The nonvanishing 
Coulomb interaction opens the possibility to achieve a 
pure spin current by tuning the orbital levels in the weak 
tunnel-coupling limit. Furthermore, we find that a pure 
spin current is obtained independently of the orbital level 
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FIG. 7. (Color online) Pumped charge and pumped spin in 
the U — ^ and the U — oo limit depending on the strength of 
the SO coupling. The illustrated functions are the maximum 
value of the pumped charge (spin) for given SO strength, aso, 
in the (e, Ae) parameter space. The coupling parameters are: 
Fl — Fr, — 7r/2, and (/)r — 7r/4. 



energies for a certain configuration of tunnel couplings, 
given by Fl Tr and 0l = -^r- 
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Appendix A: Diagrammatic rules 



We now specify the diagrammatic rules to calculate 

the diagrams of the kernels w^^'^^^^^i with n tunneling 

^ t X2X2 ^ 

lines. ^^'^^'^^'^^ Throughout the presented calculations, 
only the diagrams with one tunneling line W^*'^^ are nec- 
essary. 



1. Draw all topologically different irreducible dia- 
grams with n tunneling lines and the dot eigen- 

states Xi corresponding to Each seg- 

^ X2X2 

ment of the upper and lower contour separated by 
vertices is assigned with the corresponding eigenen- 
ergy E^(t). Each tunneling line is labeled with the 
lead A, spin <j and energie uo. 

2. Each time segment of the diagram between two 
vertices at the times tj and tj+i leads to a contri- 
bution l/{/S.Ej +zO+), where /S.Ej is the difference 
of left going energies minus right going energies. 

3. Each tunneling line contributes a factor pf{uj) or 
p(l — f{oo))^ for tunneling lines that are backward 
or forward oriented due to the Keldysh contour. 
The density of states, p, is flat and equal for the 
different leads. The Fermi function is indicated 
by f{uj) equal for the different leads as we con- 
sider no difference in the chemical potentials of the 
leads. Furthermore, each tunneling line contributes 

Ex,x' ^Ax'^aV (^fl 4 l^f) (^il ^x' l^i)' ^^^^^ I^Kf)) 
is the eigenstate in front of the vertex where the 

tunneling line starts (ends), while Xi(f)^ is the 

eigenstate of the dot behind the vertex on the 
Keldysh contour. 

4. The overall prefactor is (— 1)^+^, where b is the 
number of vertices on the lower contour line and c 
the number of crossings in the tunneling lines. 

5. Integrate over all energies of the tunneling lines and 
sum over A and a. Sum up all contributing dia- 
grams. 
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